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Abstract 

A unified canonical operator formalism for quantum stochastic differential equa- 
tions, including the quantum stochastic Liouville equation and the quantum Langevin 
equation both of the Ito and the Stratonovich types, is presented within the frame- 
work of Non-Equilibrium Thermo Field Dynamics (NETFD). It is performed by 
introducing an appropriate martingale operator in the Schrodinger and the Heisen- 
berg representations with fermionic and bosonic Brownian motions. In order to de- 
cide the double tilde conjugation rule and the thermal state conditions for fermions, 
a generalization of the system consisting of a vector field and Faddeev-Popov ghosts 
to dissipative open situations is carried out within NETFD. 
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1 Introduction 



In this paper we study time-dependent behavior of non-equilibrium quantum 
systems involving stochastic forces which can be boson or fermion type and 
are called quantum Brownian motion. Present consideration is an extension 
of previous analysis reported comprehensively by one of the authors [1] and is 
given in terms of Non-Equilibrium Thermo Field Dynamics (NETFD) [2,3,4]. 
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NETFD is a unified formalism, wliicli enables us to treat dissipative quan- 
tum systems by the method similar to usual quantum mechanics or quantum 
field theory, which accommodates the concept of the dual structure in the 
interpretation of nature, i.e. in terms of the operator algebra and the repre- 
sentation space. The representation space in NETFD is composed of a direct 
product of two Hilbert spaces: one is for non-tilde fields, and the other for 
tilde fields. Within the statistical operator (density operator) formalism there 
is entanglement between operators and statistical operator due to their non- 
commutativity. Introduction of two kinds of operators, without tilde and with 
tilde, made it possible to resolve the entanglement between relevant operators 
and the statistical operator. 

We are deriving a unified system of quantum stochastic differential equations 
(QSDEs) under the influence of quantum Brownian motion, including the 
quantum stochastic Liouville equation and the quantum Langevin equation. 
The quantum Fokkcr-Planck equation is derived by taking the random average 
of the corresponding stochastic Liouville equation. The relation between the 
Langevin equation and the stochastic Liouville equation, as well as between the 
Heisenberg equation for operators of gross variables and the quantum Fokker- 
Planck equation obtained here, is similar to the one between the Heisenberg 
equation and the Schrodinger equation in quantum mechanics and field theory. 
Our extension of analysis [1] consists of essentially three items. Two of them 
include definition of fermionic Brownian motion and treatment of fermions in 
NETFD, i.e. the tilde conjugation rule and the thermal state conditions in 
the case of fermion systems. Third item is the simultaneous consideration of 
hermitian and non-hermitian interaction Hamiltonians. 

To begin with, we first remind bricfiy some standard steps that people usually 
take in order to obtain the irreversible evolution of macroscopic systems start- 
ing from the microscopic level. At present, there are many viewpoints giving 
us tools how to describe A^-body systems out of equilibrium. At the same 
time, one usually follows one of several basic approaches: (i) the behavior of 
the systems is expressed in terms of not the total (A^-particle) distribution 
function but s-particle ones (with s being usually 1 and/or 2), (ii) the dy- 
namics of the systems is characterized by the evolution of a "coarse grained" 
phase-space distribution function or statistical operator, and (iii) the evolu- 
tion of the systems is described by the equations of motion for the dynamical 
gross variables. 

The approach (i) is intimately related to the Bogohubov method of a reduced 
description of many-particle systems [5], which is widely used for construc- 
tion of kinetic equations based on the Liouville or the Liouville- von-Neumann 
equation. Bogoliubov's hypothesis that the time dependence of higher-particle 
distribution functions enter through the one-particle distribution provides a 
fundamental importance in various schemes of truncation of the BBGKY hi- 
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erarchy. 



In the approach (ii), the most frequently used tools are projection operators 
introduced by Nakajima [6] and Zwanzig [7,8]. The basic idea underlying the 
application of their techniques to complex systems is to regard the operation 
of tracing over the environment as a formal projection in the space of the total 
system. It became especially popular in quantum optics where the so-called 
quantum master equation for reduced statistical operator of a relevant system 
now bears their names and is called the Nakajima- Zwanzig equation [9] . 

The general framework, called sub-dynamics, at the Brussels school is also 
related to the approach (ii) but the underlying concept is different from the 
one by Nakajima and Zwanzig. The main point here is the notion of the 
increase of the number of correlations within a system in time. It has been 
expounded in detail by Prigogine and coauthors, see e.g. [10,11]. 

Regarding to the approach (iii), we should mention projection operator by 
Mori [12] and the one by Kawasaki and Gunton [13]. The former is used to 
derive linear equations of motion for gross variables out of non-linear equa- 
tions. Originally, one of the intentions to introduce such an operator was to 
obtain expressions for physical (measurable) kinetic coefficients. The latter is 
an improved version of the time-dependent projection operator by Robertson 
[14]. 

Zubarev introduced the concept of non- equilibrium ensemble as a generaliza- 
tion of Green's works on the statistical mechanics of linear dissipation pro- 
cesses [15,16] and Kubo's theory of linear response of systems to mechanical 
[17] and thermodynamical [18] external perturbations. This generalization is 
known as the method of non-equilibrium statistical operator [19]. It is shown 
that this has a close relationship to the projection operator methods [20]. 

Above mentioned methods do not exhaust the entire list, but they may be the 
most generic ones. However, in this paper we do not follow them. In the case of 
presence of additional dcgrce(s) of freedom, e.g. stochastic force(s), description 
may be given also in some optional way (in a sense that consideration does 
not start from the very microscopic level) . The theory of Brownian motion is 
an example. The fundamental equation here is the Langevin equation and it 
is the stochastic differential equation for dynamical variables [21,22]. Random 
forces in Langevin equation are usually described by Gaussian white stochastic 
processes. Stochastic integral with respect to such processes is defined as a kind 
of a Riemann-Stieltjes one [23] where multiplication between the stochastic 
increment and integrand is commonly considered in the form of Ito [24] or 
Stratonovich [25] (for Ito and Stratonovich multiplications see Appendix A). 

The Langevin equation can be used to calculate various time correlation func- 
tions. Now it is radically extended to solve numerous problems in different 
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areas [26,27,28,29,30]. In particular, the theory of Brownian motion itself has 
been extended to situations where the "Brownian particle" is not a real parti- 
cle anymore, but instead some collective properties of a macroscopic system. 
Corresponding equation in the phase space or the Liouville space of statistical 
operators can be considered as a sort of stochastic differential equation too. 
In order to investigate classical stochastic systems, the stochastic Liouville 
equation was introduced first by Anderson [31] and Kubo [32,33,34]. 

There were several attempts to extend the classical theory (both the Langevin 
and the stochastic Liouville equations) for quantum cases. Study of the Langevin 
equation for quantum systems has its origin in papers by Senitzky [35,36,37], 
Schwinger [38], Haken [39,40,41,42] and Lax [43], where they investigated a 
quantum mechanical damped harmonic oscillator in connection with laser sys- 
tems. In particular, it was shown that the quantum noise, i.e. the spontaneous 
emission, can be treated in a way similar to the thermal fluctuations, and 
that the noise source has non-zero second moments proportional to a quantity 
which can be associated with a quantum analog of a diffusion coefficient. As 
it was noticed by Kubo [44] in his discussion with van Kampen, the random 
force must be an operator defined in its own Hilbert space, which does not 
happen in classical case since there is no consideration of space for the random 
force. 

Mathematical study of the quantum stochastic processes was initiated by 
Davies [45,46], Hudson [47,48,49,50,51,52], Accardi [53,54], Parthasarathy [55,56,57] 
and their co-authors. Quantum mechanical analogs of Wiener processes [47] 
and quantum Ito formula for boson systems [48.49,50,51] were defined first 
by Hudson et al. The classical Brownian motion is replaced here by the pair 
of one-parameter unitary group authomorphisms, namely by the annihilation 
and creation boson random force operators with time indices in the boson 
Fock space, named quantum Brownian motion. Fermion stochastic calculus 
were defined by Applebaum, Hudson and Parthasarathy [58,59,60,61,62,63]. 
In these papers, they developed the fermion analog of the corresponding bo- 
son theory [49,50] in which the annihilation and creation processes are fermion 
field operators in the fermion Fock space. Within the frame of this formalism, 
the Ito-Chfford integral [64,65,66,67]- fermion analog of the classical Brownian 
motion - is contained as a special case. It should be noted, however, that in 
both boson and fermion theories of quantum stochastic calculus mathemati- 
cians were debating unitary processes only. For readers' convenience, clue of 
mathematicians' theory of quantum Brownian motion and the extension with 
allowance for thermal degree of freedom are put into Appendix B. 

Contrary to expectations, attempts to extend the classical stochastic Liou- 
ville equation for quantum case were not very successful so far. In present 
work we construct our consideration using the formalism of NETFD. It is an 
alternative way to the above mentioned general methods of non-equilibrium 
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statistical mechanics in the sense that it provides us with a general structure 
of the canonical operator formalism for dissipative non-equilibrium quantum 
systems without starting from the microscopic description, and turns out to 
be especially successful in the inclusion of quantum stochastic forces. In par- 
ticular, a unified canonical operator formalism of QSDEs for boson systems 
was constructed first within NETFD [68,69,70,71,72,73,1] on the basis of the 
quantum stochastic Liouville equation. 

The paper is organized as follows. First, in section 2, we remind a brief essence 
of the formalism of NETFD by giving its technical basics and some fundamen- 
tals. In section 3 we derive the semi-free time evolution generator for systems 
in non-stationary case. The semi-free generator is bi-linear and globally gauge 
invariant. The annihilation and creation operators are introduced by means of 
a time-dependent Bogohubov transformation. We close the section by calcu- 
lating the two-point function. The generating functional method, which gives 
us the relation between the method of NETFD and the one of the Schwinger 
closed-time path, is introduced in section 4. Interaction with external fields is 
considered in section 5. Here we study two cases: hermitian and non-hermitian 
interaction hat-Hamiltonians. To make possible their simultaneous consider- 
ation we introduce an auxiliary parameter A which plays the role of a switch 
between the cases. In section 6, the general expression of the stochastic semi- 
free time evolution generator is derived for a non-stationary Gaussian white 
quantum stochastic process by means of the interaction hat-Hamiltonian with 
arbitrary A. Correlations of the random force operators are also derived gener- 
ally. With the generator, quantum stochastic Liouville equations and quantum 
stochastic Langevin equations of both Ito and Stratonovich types of the system 
are investigated in a unified manner. We conclude the section by deriving the 
equation of motion for the expectation value of an arbitrary operator of the 
relevant system. In section 7, we consider a semi-free system with a stationary 
process and check explicitly the irreversibility of such a process in terms of its 
Boltzmann entropy. In section 8 we investigate relation to the Monte Carlo 
wave-function method. Summary and open questions are put into section 9. 
Auxiliary material is put into Appendices. 



2 Basics of NETFD 



Information about the general method of NETFD can be found in many papers 
and we refer first of all to the original source [2,3,4] and the review article [1]. 
To make our paper self-contained, we include some standard steps which are 
necessary at least to fix the notations. The formalism of NETFD is constructed 
upon the following fundamental requirements. 

An arbitrary operator A in NETFD is accompanied by its tilde conjugated 
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partner A, called tilde operator, according to the rule 



(AA2)~=iii2, (2.1) 

(ci^i + 02^2)- = c^ii + (2.2) 
(i)~ = A (2.3) 

where Ci and C2 are c-numbers. It should be noted that in the present paper 
the double tilde conjugation rule (2.3) is of the same form for both bosonic 
and fermionic operators and leaves them unchanged. 

To indicate commutation or anti-commutation of two operators, say Ai and 
A2, we will use the notation [^1,^2} and call it (anti-) commutator, which 
should be understood as 

[Ai, A2} = [Ai, A2U = {Ai, A2} = A,A2 + A2A1 (2.4) 
when both operators are fermionic, or 

[A,, A2} = [Ai, A2]_ = [^1, A2] = A,A2 - A2A, (2.5) 
otherwise. ^ 

Tilde and non-tilde operators, say Ai and A2, are supposed to be mutually 
(anti-)commutative at equal time, i.e. 

[A,i2} = 0. (2.6) 

Tilde and non-tilde operators are related with each other through the thermal 
state condition (TSC) 

{9\A^^T*{9\A, (2.7) 

where (^| represents the thermal bra-vacuum; r is the complex parameter 
which takes two values: 



T— < 



1 for bosonic operators, 

(2.8) 

i for fermionic operators. 



^ When one operator is bosonic and another one is fermionic, the rule of com- 
mutation depends on the system. In this paper for such combinations we assume 
(2.5). 



6 



Derivation of the double tilde conjugation rule and TSC for fermionic opera- 
tors used in this paper is given in Appendix C. 



Within the framework of NETFD, the dynamical evolution of a system is 
described by the Schrodinger equation (here we use the system with h—1) 



^\0{t)) = -tH\0{t)), (2.9) 

where \ 0{t)) represents the thermal ket-vacuum. It can be also called the quan- 
tum master equation or the quantum Fokker-Planck equation in this paper. 
The thermal vacuums are tilde invariant, i.e., {6\'" = {6\ and \ 0{t))'" = \ 0{t)), 
and are normalized as {6\0{t)) = 1. The hat-Hamiltonian H, an infinitesimal 
time-evolution generator, satisfies the tildian condition: 



{iHy = iH. (2.10) 

The tildian hat-Hamiltonian is not necessarily hermitian operator. It has zero 
eigenvalues for the thermal bra- vacuum 



(^1^ = 0, (2.11) 
which is nothing but manifestations of conservation of probability. 
Introducing the time-evolution operator V{t) by 



!Lv{t)^-iHV{t) (2.12) 
with the initial condition y(0) = 1, we can define the Heisenberg operator 



A{f) = V-\t)AV{t) (2.13) 
satisfying the Heisenberg equation for dissipative systems 



jA{t)=i[H{t),A{t)l (2.14) 

where H{t) is the hat-Hamiltonian in the Heisenberg representation. The exis- 
tence of the Heisenberg equation of motion for coarse grained operators enables 
us to construct a canonical formalism of the dissipative quantum field theory. 
Note that with the help of TSC we have an equation of motion for a vector 

{e\A{t) 
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j^{e\A{t)=z{emtlA{t)] (2.15) 

in terms of only non-tilde operators. The expectation value of an observable 
operator A at time t is given by 



{A{t))^{e\A\0{t))^{e\A{t)\0), (2.16) 

where \0) = \ 0{t = Q)) . We define that observable operators consist only of 
non-tilde operators. ^ 



3 Semi-free hat-Hamiltonian 



Let us consider a system specified by the total hat-Hamiltonian 



Hi''' = Ht + Hi + H,,t, (3.1) 

where Ht is a semi-free hat-Hamiltonian, whereas Hi and Hi^t ^-re, respectively, 
the interaction hat-Hamiltonian within the relevant system and the one rep- 
resenting the coupling with external fields. The system itself is supposed to 
be consistent with all the requirements of NETFD given in the previous sec- 
tion. Some general remark about treatment of interaction within the relevant 
system is given in section 4. Explicit treatment of interaction with external 
fields is given in section 5. Here we concentrate on derivation and study of 
properties of the semi-free hat-Hamiltonian, i.e. renormalized unperturbed 
hat- Hamilt onian . 



3.1 Derivation of the semi-free hat-Hamiltonian 



The semi-free hat-Hamiltonian is bilinear in operators a, a\ a and a\ and is 
invariant under the phase transformation a — > ae"^: 



H^ = hi{t)a^a + h2{t)a)a + h^{t)aa + h^{t)a^d^ + hoit), (3.2) 

where hj{t) are time-dependent complex c-number functions. Operators a, a\ 
a and satisfy the canonical (anti-) commutation relations 

^ We can include tilde operators in addition to non-tilde ones in the definition of 
observable. However, inclusion of tilde operators may give us a set of different but 
equivalent definitions for one observable operator. 



8 



(3.3) 



where we use a = 1 for bosonic systems and a = —1 for fermionic ones. 
According to (2.6), tilde and non-tilde operators are mutually (anti-)com- 
mutative. In the following account, a subscript k for specifying a momentum 
and/or other degrees of freedom will be dropped unless it is necessary. Number 
of unknown functions hj(t) can be reduced by the use of (2.11) and tildian 
(2.10) for the semi-free hat-Hamiltonian. It results in 



where 



(3.4) 



nt = ci(t) (a^a + a^a) -r[2ci(t) + C2{t)]o)a^ 
+ aTC2{t)aa + a [2ci(t) + C2{t)] , 



with 



(3.5) 
(3.6) 



u{t)=^thi{t), 
ci(t)^ Qmhi{t), 



(3.7) 
(3.8) 



C2{t) = < 



^mhs{t) for bosonic systems, 
^ths{t) for fermionic systems. 



(3.9) 



Let us introduce operators a{t) and a^{t) in the interaction representation 
defined by 



a{t) = V-\t)aV{t), a^{t) = V-\t)a^V{t), 
where 



(3.10) 



±V{t) = -zHtV{t), 



(3.11) 



with the initial condition ^(0) = 1. They satisfy the equal-time (anti-)commutation 
relations 



[a(t),att(t)]_, = l. 



(3.12) 
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The Heisenberg equation (2.14) for a{t) and a^{t) with 



H{t) = V-\t)HtV{t) 
are explicitly given by 



(3.13) 



da{t) 

dt 
da^{t) 
dt 



= [ci{t) - iuit)] a{t) - T [2ci(0 + C2{t)] att(0, 
= [iu;{t) - ci{t)] att(i) - rc2{t)a{t). 



(3.14) 
(3.15) 



In these formulae we used a symbol jj instead of usual dagger because the 
semi-free hat-Hamiltonian is not necessarily hermitian. 

Since the semi-free hat-Hamiltonian Hf satisfies (2.11), we have TSC for the 
bra- vacuum at time t 



\d^{t)^r*{e\a{t). 



(3.16) 



By making use of the Heisenberg equations (3.14) and (3.15), and of TSC 
(3.16), one obtains the equation of motion for a vector {9\a^{t)a{t) in the 
form 



^{e\a^{t)a{t) = -2K{t){e\a'^{t)a{t)+iJ:<{t){e\, 
dt 

where K.{t) and iS<(t) are defined by 

/t(t)=Ci(t) +C2(t), 

iS<(t) = -a[2ci(t) + C2(t)]. 



(3.17) 



(3.18) 
(3.19) 



Substituting (3.18) and (3.19) into (3.14) and (3.15) one gets equations of 
motion for operators a{t) and a'^{t) in the form 



da{t) 

dt 
da^{t) 

dt 



= - [iuj{t) + K{t)] a{t) - (TiL<{t) a{t) - Ta^{t) 



[iu{t) + K{t)\ a^{t) + (7iE<(i) a^{t) - Ta{t) - 2TK{t)a{t)i?,. 21) 



(3.20) 



Applying the thermal ket-vacuum \0) at the initial time to (3.17), we obtain 
the equation of motion for the one-particle distribution function 
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n{t)^{9\a^{t)a{t)\0) 
in the form 



(3.22) 



^n(t) = -2K(t)n(t) + iE<(t). (3.23) 
at 

Equation (3.23) can be identified as the generahzed Boltzmann equation of the 
system. The function iE^ (t) is given when the interaction hat-Hamiltonian Hi 
is defined. 

The initial ket- vacuum 1 0) is specified by TSC 



d\0) = Tfa^\0) (3.24) 

with / e R. The initial value for the one-particle distribution function n — 
n{t — 0) is determined by /. Since 



n^n* = {9\a''a\0)^ = {9\a^a\0) = r*{9\aa\0) = \r\^f{9\aa^0) 
= /[l + H, (3.25) 

we have 



n = f[l-af]-\ (3.26) 

In the first equahty of (3.25) we used the fact that n is a real number; in the 
third equahty we used the tilde invariance of the thermal vacuums (^| and 
\0); finally, in the fourth and fifth equalities we used TSCs (2.7) and (3.24), 
respectively. 

Solving the Heisenberg equations for a{t), a''^{t) and their tilde conjugates, 
and using TSC at initial time (3.24), we find TSC for the ket- vacuum at time 
t 



m\0)^ ^l''^^^^^^ aHt)\0), (3.27) 
where n{t) satisfies the Boltzmann equation (3.23). 

Substituting (3.18), (3.19) and (3.23) into (3.6) one gets the most general form 
of ilt in the interaction representation: 



flt^- {K{t) [1 + 2an{t)] + ah{t)} (a^a + a^a) 
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+(7r {2K{t) [1 + an{t)] + an{t)} ad 

+aT {2K{t)n{t) + h{t)} a^a^ - {2K{t)n{t) + h{t)} . (3.28) 
Here, we used the abbreviation n{t) = dn{t)/dt. 
By introducing thermal doublet notations 

^ (a\ -Toj , a'' = coUon (a, ra^) , (3.29) 
canonical (anti-)commutation relations are written as 

(3.30) 

The resulting semi-free hat-Hamiltonian (3.4) can be presented in a compact 
form as 

Ht^u;{t)d^a^' + id^A{t)^''a'' + a[u;{t) + iK{t)], (3.31) 
where matrix A{tY'^ has the following structure: 



^ ^ I [2n(i) +a]- h{t), 2K{t)n{t) + h{t) 

-2K{t) [n{t) + 0-] - h{t), K{t) [2n{t) + cr] + n{t) 



3.2 Annihilation and creation operators 

Let us introduce annihilation and creation operators by 

ry[t) = [1 + (Tn{t)\ a{t) - (TTn{t)a^{t), (3.33) 

')^{t)^a^{t)-aTa{t). (3.34) 

From TSCs (3.16) and (3.27) at time t, we see that they annihilate the vacu- 
ums: 

(^|7^(t)=0, i{t)\0) = 0, (3.35) 

(^|7^(t) = 0, mO) = ^. (3.36) 

With the thermal doublet notations 
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a{tY = (att(i), -Ta{t)) , a{ty = coUon (^a{t) , ra^ (t)) , (3.37) 
7W^= (7^(t), -r7(t)) , = coUon (7(t), r7^(t)) , (3.38) 

(3.33), (3.34) and their tilde conjugates can be written as 

7(i)" = a{ty[B-\t)]'"' , 7(i)^ = BitY^aity, (3.39) 
where B{t)^^ is a matrix of the time-dependent Bogoliubov transformation: 



B{t) 



1 + an{t) -an{t) ^ 



-1 1 



(3.40) 



This transformation is the canonical one since it leaves the canonical (anti-)- 
commutation relations unchanged: 

[j{tr,j{tru^s^\ (3.41) 

The equation of motion for the thermal doublet ^{tY is derived as 

j^^itr = Huj{t)6'^'' - K{t)Tn 7(^)^ (3.42) 

where 

A 



1 
-1 



(3.43) 



The solution of (3.42) then is obtained in the form 



^{ty = |/ i-Mt')^'"" - <t')Tr] I 7(0)^ (3.44) 



3.3 Schrddinger representation 



Annihilation and creation operators in the Schrodinger representation are in- 
troduced by the relations 
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^^tr^v-\t)^:^v{t), ^{tr^v-\t)^^v{t), (3.45) 

with V{t) being specified by (3.11) and the thermal doublet notations 

It = (7^, -TJt) , It = coUon {yt, t7^) ■ (3.46) 
Using (3.29), one can write 

-f!^ = a^[B-\t)r, ^^ = B{tY''a\ (3.47) 
where matrix B{tY^ is given by (3.40). 

We see that the annihilation and creation operators in the Schrodinger repre- 
sentation annihilate the vacuums at time t: 

(^|7^ = 0, 7,|0(t))=0, (3.48) 
(^|7^ = 0, ^^\0{t))=Q. (3.49) 

Note that creation operator 7"^ and its tilde conjugated partner 7"^ do not 
depend on time. It is consistent with the fact that the vacuum (^| does not 
depend on time due to the property {0\Ht = 0. 

The terms (3.5) and (3.28) of the semi-free hat-Hamiltonian Ht, (3.4), now 
read 

^s,t = ^i) (7^-7^7*), (3.50) 
nt = -Kit) (7*7i + 7*7^) + (7Tn(t)7*7^ (3.51) 

in the normal ordering with respect to the annihilation and creation operators 
in the Schrodinger representation. When the system is semi-free, putting (3.50) 
and (3.51) into (3.4) and substituting Ht for H into the Schrodinger equation 
(2.9), one has 



^^\0{t))^aTn{tY^^^\0{t)). (3.52) 
It is solved to give 



\0{t))^ex^{-{e\^tlt\0)l'^l^}\0). (3.53) 
Here we introduced a kind of the order parameter 
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{e\%lt\0)^aT[n{Q)-n{t)] 



(3.54) 



which gives a measure of difTerence of the system from the initial state. Prom 
(3.53) we see that the evolution of the ket- vacuum is realized by a condensation 
of tilde and non-tilde particle pairs into initial kct- vacuum. The ket-vacuum 
itself is the functional of the one-particle distribution function n{t). The de- 
pendence of the thermal ket-vacuum on n{t) is given by 



'0{t))^aT-i^^^\0{t)). (3.55) 



5n{t) 

Then the Schrodinger equation can be written in an alternative way: 



This shows that the vacuum \0{t)) is migrating in the super-representation 
space spanned by the onc-particle distribution function {nk(t)} with the ve- 
locity {h\f_{t)} as a conserved quantity [74,75]. 



3.4 Two-point function of the semi-free field 

A time-ordered two-point function G{t,t'Y^ (propagator), defined by 

G{t, t'y = -i{e\T[a{tYa{t'Y\\0), (3.57) 
is given by 

G(t, t'^" = [B-\t)Y^g{t, fY^Bity, (3.58) 

with 



g{t,t'Y^ ^ -z{e\T[j{tYj{tr]\o) 

where non-zero matrix elements are 
G^{t,t')^-i{9\TMth^{t')]\0) 



^G^(t,t')0 ^ 
G'^{t,t')J 



(3.59) 



-i9{t - t') exp dt" [-iuj{t") - n{t")] I , (3.60) 



15 



G^{t,t')^ia{0\T[^^mt')]\O) 

t' 



■ ie{t' -t)exp < 



Jdt" [iuj{t") - K{t") 



> . 



(3.61) 



Here, T is the time ordering operator. 



4 Generating functional method 



Let us define the generating functional for the semi-free field by 



i / dt S{t) 



Z[K,K] = {9\Tcxp< 



where the source function S{t) reads 



(4.1) 



S{t) = K{tYa{tY + d{tYK{tY = K^{tYi{tY + i{tYK^{tY. (4-2) 



with 



K{tY = {K{t)\ -rK{t)) , K{tY = coUon [K{t),rK{t)*) , (4.3) 

and similar notations for K^itY and K^{tY- The K^s are related by the Bo- 
goliubov transformation 



K^{tY = KitYlB-^it)^, K^{tY = B{tY"K{tY. 



(4.4) 



Matrix B{tY^ here is the one given by (3.40). External fictitious fields K{tYi 
K(tY are c-numbers or Grassmann numbers corresponding to cr = 1 or cr = 
— 1, and satisfy 



[K{tY,K{tYU-^- 



(4.5) 



Operators a{tY-i ciiiY i li'^Y "^{tY are those in the interaction represen- 
tation introduced in section 3. 

Taking the functional derivative of the generating functional (4.1), one has 



5\TiZ[K,k] = ~ijdt [5K,{tY{imK + mY)JK,{tYl (4.6) 
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where (7(i)^)^ and (7(i)'^)^ are defined by 



1 



lnZ[K, k] 



^r^ ^MTbitr exp{-t / dt'S{t')}]\0), 



5 



In Z[K, k] 



1 }' 
^{e\mtreM-ij dt'smio). 



The equation of motion for (7(i)^)^ is obtained in the form 



(4.7) 



(4.8) 



With the boundary conditions 

(7(0)^=^), = (7(0)), = 0, 
(7(tr=').=r(7"(t)).=0, 

(7(0)-='), = -r(7(0)),=0, 
it can be solved as 



(4.9) 



(4.10) 
(4.11) 
(4.12) 
(4.13) 



{imK=Jdt'g{t,tTK,{tr, 



(4.14) 



where Q{t.t'Y^ is given by (3.59). The boundary conditions in (4.10) to (4.13) 
are derived by TSCs (3.35) and (3.36). 

Substituting (4.14) into (4.6), one finally obtains [76] 



Z[K, X]=exp< 



: exp < 



t t 

-i Jdt J dt' K^{tYg{t,tyK^{t'y 



t t 



-ijdtj dt' k{tYG{t,tyK{t'y 







> . 



(4.15) 
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This expression has been derived first by Schwinger for a boson system within 
the closed-time path method [38] . Derivation of this result shown in the present 
section reveals the relation between the quantum operator formalism of dissi- 
pative fields (reahzed for the first time within NETFD) and their path integral 
formalism [38]. 

The effect of the interaction Hi = Hi{a^,a'^) within the system, which in- 
duces the dynamical correlations, can be taken into account by the generating 
functional 



>Z[K,K]. (4.16) 




Note that Hi should satisfy {e\Hi = 0. 

5 Interaction with external fields 

5.1 Hermitian interaction hat-Hamiltonian 



The simplest hat-Hamiltonian representing an interaction with an external 
field may be given by 



H[^H[-H[, (5.1) 
with a hermitian interaction Hamiltonian 

H'^^i(a\-h\a), (5.2) 

where 6^, h\ and their tilde conjugates are operators of the external system 
and are assumed to be (anti-) commutative with operators a, and their 
tilde conjugates of the relevant system. The subscript t indicates that these 
operators may depend on time. Note that the hat-Hamiltonian (5.1) is tildian, 
i.e. 

{iH^y^iH',. (5.3) 

The tilde and non-tilde operators of the external system are related with each 
other by 
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{\bl-r*{\h, (5.4) 

where (| is the bra- vacuum for the external system. Applying the bra- vacuum 
(6*1 for the relevant system on (5.1), one has 



{e\Hi = -z{e\((3fa + aTpfa^). (5.5) 
Here we introduced a new operator 

P! = bl-Tbt, (5.6) 
which annihilates the bra- vacuum (| for the external system: 

(1/3/ = 0. (5.7) 

As it is seen from (5.6), the subscript t of the new operator has been 
inherited from the original operators of the external system. By applying the 
bra- vacuum (| on in addition to (^|, we observe that 

mm = 0, (5.8) 
where the bra- vacuum of a total system is introduced by 

m^{\-{e\. (5.9) 



The dynamics of the system is described by the Schrodinger equation for the 
ket- vacuum \0{t))) of the whole system: 



where Hi^t in H^°^ is replaced by H^. Conservation of the probabihty is guar- 
anteed by {{9\H^°^ — for the total system, i.e. the relevant system and the 
external system. 



5.2 Non-Hermitian interaction hat-Hamiltonian 

Let us consider if we can have an interaction hat-Hamiltonian which satisfies 
the conservation of probability within the relevant system. This feature is 
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consistent with the one we have in the case of stochastic differential equations 
for classical systems. 

We assume that the interaction hat-Hamiltonian is globally gauge invariant 
and bilinear: 

H" = i ^hia^bt + h2a^bl + hsabt + h4dbl 

+ h^a^bt + hgd^bl + hrobt + hgabl^ , (5-11) 

where quantities hj {j — 1, ■ ■ ■ ,8) are time- independent complex c-numbers. 
The tildian 

{iH^y^iH^' (5.12) 
gives us 

hi = /i5, = h^, hi = hj, hi = kg. (5.13) 
By applying (^| from the left to the Schrodinger equation 



with Hj^t in H^°^ being replaced by H^', we see that the requirement of the 
conservation of probability within the relevant system leads to 



{e\H't' = 0. (5.15) 

Ht in is the semi-free hat-Hamiltonian of the relevant system satisfying 
(2.11). From (5.15) we obtain 



hi + arks = 0, hj + arh^ = 0, (5.16) 
/t2 + crrh^ = 0, hs + arh^ = 0, (5-17) 

which are solved as 



hs — —rhi, h-j — —arhl, (5.18) 
/i4 = — t/i2, ^8 — — cr/i2- (5.19) 

Then the structure of H" can be expressed in terms of only hi, /i2 and their 
complex conjugates as 
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H^'^i{a^Pt + a^Pt}: (5.20) 
where we introduced new operators 

= a^ — ra, (5-21) 
Pt = hibt + h2bl (5.22) 

and their tilde conjugates. Note, that the creation operator annihilates the 
bra- vacuum (^|: 

{e\a^^O. (5.23) 
In order to investigate parameters hi and /i2 we consider the moments 

(AA) = {h + arh) {rhliblh) + h;{bM)} , (5.24) 
(AA) = {hi + Th;) {arhiiblbt) + h2{btbl)} , (5.25) 

where we are using the symbol (• • •) = (| • • • |^) without specifying the dy- 
namics which determines the ket-vacuum |t) of the external system. For the 
present purpose, the details of its dynamics are not required. Here we assume, 
however, that the external ket-vacuum may evolve in time. The further use 
of the property of the (anti-)commutativity, i.e. (AA) = cr(AA)) gives the 
necessary two relations to define hi and h2: 

(rhi + h2)hl^(rhl + ah*2)hi, (5.26) 
alrhi + h2)hl = (r/i* + ahi)h2, (5.27) 

which reduce to 

hlh2 = ahihl. (5.28) 
We can express hi and /i2 as 

hi^lie'^\ h2^ve'^\ (5.29) 

where /x, e R, namely ii = \hi\, v = \h2\. From the requirement (5.28), one 
has 92 — Oi lor a — 1 and 62 — Oi — 7r/2 for a — —1. Substituting (5.29) into 
(5.22) and putting the phase factor e*^i into bt and foj, we have 

Pt^f^bt + ariybl (5.30) 
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Thus, the vector {\(3t is calculated as 



{\l3t^{\{lxht + aTvh\)^{li + av){\ht. 



(5.31) 



The further requirement that the norm of {\(3t should be equal to that of 
i.e. \\{\(3t\\ = IKI&tll, leads one to the relation 



5.3 Relation between the two interaction hat-Hamiltonians 

Note that the hcrmitian interaction hat-Hamiltonian H[ and the non-hermitian 
one Hi are related to each other by 



With an auxiliary parameter < A < 1, it is possible to make a simul- 
taneous consideration of both hermitian and non-hermitian interaction hat- 
Hamiltonians by introducing 



It is easy to see that this expression is reduced to H[ in the case A = 1 and 
to H'{ in the case A = 0, respectively. The dynamics of the system is now 
described by the Shrodinger equation (5.10) or (5.14) with Hi^t in Hl°^ being 
given by (5.34). 



6 Quantum stochastic differential equations 

6.1 Quantum stochastic Liouville equations 
6.1.1 ltd type 

Let us derive the general form of the semi- free hat-Hamiltonian li.F,tdt for a 
stochastic Liouville equation of the Ito type 



fj, + au — 1. 



(5.32) 




(5.33) 



arVa^)}. (5.34) 



d\OF{t)) ^ -inF,tdt\OF{t)) 



(6.1) 
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where a subscript F is added to indicate that we are considering a system 
under the influence of a random force. We assume that the hat-Hamiltonian 
'HF,tdt for the stochastic semi-free field is bilinear in a, a\ dF^, dF} and their 
tilde conjugates, and that it is invariant under the phase transformation a — > 
ae'*^ and dFt — > dFfB^'^. Here, a, a' and their tilde conjugates are operators of 
a relevant system satisfying the canonical (anti-)commutation relation 



kat]_, = l, (6.2) 

whereas dFt, dF} and their tilde conjugates are random force operators. The 
tilde and non-tilde operators are related with each other by the TSC 



{e\~a) = T*{e\a, (6.3) 
{\dF}^T*{\dF,, (6.4) 

where (^| and (| are, respectively, the thermal bra- vacuum of the relevant 
system and of the random force. 

From the investigation in section 5, we can propose that the required form of 
the hat-Hamiltonian should be 



'HF,tdt = Hfdt + dMt, (6.5) 

where Ht is specified by (3.4) with fit having the same structure as (3.28) or 
(3.51). For later convenience, we rewrite iJ* as 



ilt^il^^il^ (6.6) 

with 

-^R = — '^(i) {tt^a + a^aj , (6.7) 

nY)^(TT{2n{t)[n{t)+ri\+n{t)}a^a^, (6.8) 

where we introduced 

a = ^a + aT'qd\ ^ + a'q=l, (6.9) 
which forms a canonical set with defined by (5.21), i.e. 

[a, = 1. (6.10) 
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The one-particle distribution function n{t) is defined by 



it) = ({e\a^a\OF{t))), (6.11) 



and satisfies the Boltzmann equation (D.ll) (sec Appendix D). Here, (■ • •) 
means to take the random average, i.e., the vacuum expectation value with 
respect to the thermal bra- and ket- vacuums of random force: (■■•) = (| ■■■ |)- 
Terms IJ^ and TJd are, respectively, the relaxational and diffusive parts of the 
damping operator IJt. 

The martingale dMt is the term containing operators representing quantum 
Brownian motion and satisfies 



{dMt) = Q. (6.12) 
Associating dFt and dF^ with bfdt and bldt in (5.34), respectively, we have 

dMt = i {a'^dWt + a'^dWt} 

-iX {dW^ina + ctW) + dW^^i^d + arVa^)} , (6.13) 

where we introduced new operators 

dWt = i^dFt + (TTvdF} , (6.14) 
dW^^dF} -TdPf (6.15) 

Note that dW^ and dWt annihilate the bra- vacuum for random force (|: 

(|ciW/ = 0. (6.16) 

In the Ito multiplication, the random force operators dWt, dWt and their tilde 
conjugates do not correlate with quantities at time e.g. \OF{t)): 



dMt\OF{t))) ^0. (6.17) 

Thus, taking the random average of the stochastic Liouville equation (6.1), 
we arrive at the Fokker-Planck equation 



^^\0(t))^-iHt\0(t)), (6.18) 
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where \ 0{t)) = {lOpit))). 



The formal solution of (6.1) can be written as 



OF{t))^VF{t)\OFm, 



(6.19) 



where the time-evolution generator is defined through 



dVpit) = -iHF,tdtVF{t) 



(6.20) 



with the initial condition Vf(0) = 1. 

6.1.2 Fluctuation- dissipation theorem of the second kind 

By making use of the relation between the Ito and the Stratonovich stochas- 
tic multiplications (see Appendix A), we can rewrite the Ito type stochastic 
Liouville equation into the Stratonovich type as follows. Relation (A. 6) makes 
the term containing random force operators in the r.h.s of (6.1) be 



where the symbol o has been introduced to indicate the Stratonovich stochastic 
multiplication (see Appendix A). Substituting (6.1) into the last term for 

d\OF{f^) and neglecting terms of the higher order than d.t. we arrive at the 
quantum stochastic Liouville equation of the Stratonovich type 



In order to obtain expression (6.24) we used the generalized fluctuation- 
dissipation theorem of the second kind, which can be written as 




(6.21) 



d\OF{t))^-iHF,tdto\OF{t)), 



(6.22) 



with 



HF,tdt = Hs,tdt + iiltdt + dMt + -dMtdMt 
= Hs,tdt + i{l- X)nKdt + dMf 



(6.23) 

(6.24) 




(6.25) 
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(refer to Appendix D for derivation). Note that in the Stratonovich multiph- 
cation random force operators dWt, dW^ and their tilde conjugates correlate 
with quantities at time t, i.e., 



(dMto\OF{t))) ^Q. (6.26) 

The formal solution of (6.22) has the form (6.19), where the stochastic time- 
evolution generator VF{t) is defined through 



dVpit) = -iHF,tdt o VF{t) (6.27) 
with the initial condition Vf(0) = 1. 

6.1.3 Correlations of the random force operators 

Operators dWt, dW^ and their tilde conjugates are of the quantum stochastic 
Wiener process satisfying (for derivation of the results see Appendix D) 

(dWt) = (dWt) = 0, (6.28) 

(dWtdWs) = (dWtdWs) = 0, (6.29) 

{dWtdWs)^a{dW,dWt) 

^T{2K{t)[n{t) + i']+h{t)}S{t - s)dt ds, (6.30) 

and 

{dW,^) = {dW,^)=0, (6.31) 

{dWfdW^) = {dW^^dWf) = 0, (6.32) 

{dWldWs) = {dWldWs) = 0, (6.33) 

{dWtdW^) = {dWtdWl) = 2n{t)5{t - s)dt ds. (6.34) 

Due to the argument of Appendix D we also have ^ — and rj — u, which 
leads to 

a — jjia + arudK (6.35) 

In the following, we will use this definition in both 77^ and dMf. Especially, 
the latter becomes 

dMt = i {a^dWt + a'^dWt} - iX {dW/a + dW^a} . (6.36) 
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It is important to note here that the martingale dMt is introduced in the 
normal ordering with respect to all operators a^, a, dWt^, dWt and their tilde 
conjugates. 

Within the weak relations, the correlations (6.29), (6.30) and (6.32) to (6.34) 
reduce, respectively, to 



dWtdWs = dWtdWs = 0, (6.37) 

dWtdWs = adWsdWt = r {2K{t)[n{t) + u] + h{t)} 5(t - s)ds dt, (6.38) 

and 

dW^^dW^ = dW^dWf = 0, (6.39) 

dWfdWs = dWfdWs = 0, (6.40) 

dWtdW^ = dWtdW^ = 2K{t)6{t - s)ds dt. (6.41) 

6.2 Stochastic semi-free operators 

The stochastic semi-free operators are defined by 

A{t) = Vp\t)AVF{t), (6.42) 
whereas the random force operators in the Heisenberg representation by 

Wit) = Vf\t)WtVFit), W^it) = Vf\t)WlVFit), (6.43) 

and their tilde conjugates. We also use the convenient operators introduced 
by 

aW{t) = Vi\t)dWtVF{t), dW^{t) = Vf\t)dWt^VF{t), (6.44) 
and their tilde conjugates. Here, 

dVf^ (t) = iVf\t)np^^dt, (6.45) 
with Vf^{0) — 1. Hp^tdt is specified by 

Hp tdt = Hp^tdt + idMtdMt = Hs,tdt + itfdt + dMt (6.46) 



27 



with 



i7r = (l-2A)i7R-i7D. (6.47) 
In particular cases when A represents a or we have 

a(t) = Vp\t)aVF{t), att(t) = Vf\t)a^VF{t). (6.48) 

Since the stochastic tildian hat-Hamiltonian Tip^tdt is not necessarily hermi- 
tian, we introduced the symbol ff in order to distinguish it from the hermite 
conjugation f . It is assumed that, at initial time t = 0, the relevant system 
starts to contact with the irrelevant system representing the stochastic process 
described by the random force operators dFt, dP} and their tilde conjugates. 
Within the formalism, the random force operators dFt and dP} are assumed 
to (anti-) commute with any relevant system operator A in the Schrodinger 
representation, i.e. 

[Ac^Fi} = 0, [A,dF}] = Q. (6.49) 

The semi- free operators a(i), a^{t) and their tilde conjugates keep the equal- 
time canonical (anti-)commutation relations 

[a(i),att(i)]_, = l, (6.50) 
and satisfy TSC 

{{e\~a^{t)=T*{{e\a{t). (6.51) 

Calculating the time derivatives of Heisenberg operators of the quantum Brow- 
nian motion (6.43) within the Ito calculus (A. 15), and taking into account 
(6.20), (6.45) with the characteristics of the Ito multiphcation 

[dWl, Wt]-. = [dWt, W^*]-. = 0, (6.52) 
[dWt, Mt)] = [dW,^, Mt)] = 0, (6.53) 

and their tilde conjugates, one has 

dW(t)^dWt - iVf\t)[dWt, dMt]Vp(t), (6.54) 
dW^(t)^dWl - iVi\t)[dWt^, dMtlVpit), (6.55) 
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and their tilde conjugates. With the help of (6.37) to (6.41), the expressions 
(6.54) and (6.55) reduce, respectively, to 



dW{t) = dWt - 2XK{t)a{t)dt, (6.56) 
dW^{t)^dWt^ - 2K{t)a^{t)dt, (6.57) 

(and their tilde conjugates), while (6.53) gives 

dW{t) = dWt, m^{t) = dWf, (6.58) 
(and their tilde conjugates). 

We see that the martingale operator dM{t) = Vf^{t)dMtVF{t) being written 
in terms of the Heisenberg operators reads 

dM{t) = i{a^{t)dW{t) + &^{t)dW{t)} 

- i\ [dW^{t)a{t) + dW^{t)a{t)} 
^i{a^{t)dW{t) + a^{t)dW{t)} 

- iX {dW^{t)a{t) + dW^{t)a{t)} 

= dM{t), (6.59) 

and keeps the property 

{dM(t)) ^ {dM(t))^0 (6.60) 

for arbitrary A. The beautiful relation (6.59) is manifestations of the normal 
ordered definition (6.36). Note that the increments in the martingale (6.36) 
are introduced just through the random force operators dWt', dWt and their 
tilde conjugates. Therefore, dM{t) is different from the operator calculated by 
d (i {a^{t)W{t) + t.c.} - iX \w^{t)a{t) + t.c.}). Here, t.c. indicates the tilde 
conjugation. 



6.3 Quantum Langevin equations 



6.3.1 ltd type 

Substituting dVpit), (6.20), and dVp^it), (6.45), into the time derivative of 
the dynamical quantity A(t) within the Ito calculus (A. 15), we obtain the 
quantum Langevin equation of the Ito type in the form 
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dA{t)^i[nFit)dt,A{t)]-dM{t)[dM{t),A{t)] (6.61) 
^i[Hs{t),Ait)]dt 
+K{t) (a^ (t) [a{t) ,A{t)} + a'^ (t) [a{t) ,A{t)} 

+ (2A - l){[A{t), a^{t)}a{t) + [A{t), a*(t)}a(t))) dt 

+T{2K{t)[n{t) + u] +n{t))[a^{t), [a^{t),A{t)}}dt 

+ [A{t),a^{t)}dW{t) + [A{t),a^{t)}dW{t) 

+A [dW^{t)[a{t), A{t)} + dW*{t)[a{t), A{t)}) (6.62) 

^i[Hs{t),A{t)]dt 

-K{t) ([A{t),a^{t)}a{t) + [A{t),a^{t)}a{t) 

+ {2X-l){a^{t)[a{t),A{t)} + a^{t)[a{t),A{t)}))dt 
+T{2K{t)[n{t) + u] + n{t))[a^{t), [a^{t), A{t)}}dt 
+[A{t), a^{t)}dWt + [A{t), a^{t)}dWt 

+X (^dWf[a{t),A{t)} + dWf[a{t),A{t)}) , (6.63) 

where (6.37) to (6.41) for multiplications among the random force operators 
are employed. To derive (6.63) from (6.62), we used (6.56) and (6.57). Note 
that the Langevin equations (6.62) and (6.63) written, respectively, by means 
of the quantum Brownian motion in the Heisenberg representation and by 
means of that in the Schrodinger representation may be related with the "out" 
and "in" fields introduced by Gardiner et al. [77,78]. 

With the help of (6.62) one can verify that the calculus rule for the product 
of arbitrary relevant stochastic operators, say A{t) and B{t), satisfies the Ito 
calculus (A. 15). This proves that QSDE (6.62) is of the Ito type indeed. Fur- 
thermore, since (6.62) is the time-evolution equation for any relevant stochastic 
operator A{t), it is Ito's formula for quantum systems. 



6.3.2 Stratonovich type 

The quantum Langevin equation of the Stratonovich type can be derived sim- 
ilarly if one starts from the expression for a dynamical quantity A{t), (6.42), 
and considers its derivative in the Stratonovich calculus (A. 16) with the help 
of 



dVf\t) = iVf\t) o Hp^tdt. (6.64) 

Substituting (6.27) for dVplt) and (6.64) for dVf^{t) into dA{t), we have for 
the stochastic Heisenberg equation of the Stratonovich type 
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dA{t)^i[HF{t)dt ° A{t)] (6.65) 
^i[Hs{t),Ait)]dt 
+K{t) (a^ (t) [a{t) ,A{t)} + a'^ (t) [a{t) ,A{t)} 

+ (2A - l){[A{t), a^{t)}a{t) + [A{t), a*(t)}a(t))) dt 

+ [A{t), a^{t)} o dW{t) + [A{t), a^{t)} o dW{t) 

+A [dW^{t) o A{t)} + o A{t)}) (6.66) 

= z[^s(t),^(t)]t^t 

-«;(i) ([^(i),a*(t)}a(t) + [A{t),a^{t)}a{t) 

+ (2X-l)(a^(t)[a(t),A(t)} + a^(t)[a(t),A(t)}))dt 

+[A{t), a^{t)} o dWt + [A{t), a^{t)} o dWt 

+X [dWf o [a{t),A{t)} + dWl o [a{t),A{t)}) . (6.67) 

Here, we defined 

[X(t) ° Y(t)] = X(t) o Y{t) - Y(t) o X(t) (6.68) 
for arbitrary operators X{t) and Y{t), and 

HF{t)dt = Vi\t) o Hp^tdt o VF{t). (6.69) 
Note tfiat 

Vf\t) o dMt o Vpit) = Vp\t)dMtVF{t) 

+^Vi\t)dMtdVF{t) + ^dVi\t)dMtVF{t) 
^dM{t). (6.70) 

Using expression (6.66), one can readily verify that the calculus rule for the 
product of arbitrary relevant system operators, say A{t) and -B(t), satisfies 
the Stratonovich type calculus (A. 16). This fact proves that QSDE (6.66) 
is indeed of the Stratonovich type, and provides us with the reason why the 
stochastic Heisenberg equation (6.65) has the same structure as the one (2.14) 
for non-stochastic operators. 

The quantum Langevin equation of the Stratonovich type can be also derived 
from that of the Ito type by making use the connection formulae (A. 13) and 
(A.14). When dY{t) is dW(t) or dW{t), and X(t) is constituted by the relevant 
operator, say A{t), satisfying the quantum Langevin equation of the Ito type, 
the connection formula (A. 13) reduces, respectively, to 
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Similarly, when dX(t) is dW^{t) or dW^{t), and Y{t) is A{t), the connection 
formula (A. 14) reduces, respectively, to 



dW^{t)-A{t) = dW*{t)oA{t), 
dW^{t) ■ A{t) = dW^{t) o A{t). 



(6.73) 
(6.74) 



Using these relations in (6.62), the quantum Langevin equation of the Stratonovich 
type is obtained in the form (6.66). 

Substituting a and for A as an example, we see that both (6.61) and (6.65) 
result in 



which are written in terms of the original operators as 

da{t) = - [iu;{t) - \K{t)] a{t)dt 

- (1 - X)K{t)[{ii - aiy)a{t) + 2aTud^{t)]dt 
+ dW{t) - XarudW^it) 

= - [iu;{t) + XK{t)] a{t)dt 

- (1 - X)K{t)[{tjL - (Tv)a{t) + 2(TTva^{t)\dt 

+ dWt- XarudWt^, (6.77) 
datt(i) = - [iu;{t) - XK{t)] a^{t)dt 

+ (1 - X)n{t)[{n - (Tv)a^{t) - 2(TTna{t)\dt 
+ aTdW{t) + XndW^{t) 
= - [iuj{t) + Xn{t)\ a^{t)dt 

+ (1 - X)n{t)[{n - av)a^{t) - 2aT fjLa{t)]dt 

+ ardWt + XiidW^. (6.78) 

These equations are the same in both Ito and Stratonovich multiplications as 
they should be with the martingale (6.36). 



da{t) = -iuj{t)a{t)dt - (1 - 2X)K{t)a{t)dt + dW{t) 
= - [iu;{t) + K{t)] a(t)dt + dWt, 
da^{t) = -iu{t)a^{t)dt + K{t)a*{t)dt + XdW^{t) 
= - [iu;{t) - (1 - 2X)K{t)] a^{t)dt + XdW^^, 



(6.75) 



(6.76) 
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6-4 Averaged equation of motion 



Applying the total bra- vacuum ((^| to the Ito type quantum Langevin equa- 
tion (6.61) to (6.63), one can derive the stochastic equation of motion of the 
Ito type for the bra- vector state {{9\A{t) in the form 

d{{9\A{t)^i{{e\[Hs{t),A{t)]dt 

+(2A - l)K{t){{9\ (A{t)[a^{t)a{t) + a^{t)a{t)]) dt 
+ar{2K{t)[n{t) + u] + h{t)){{e\A{t)a^{t)a^{t)dt 

+ {{9\ [A{t)[Q.^{t)dW{t) + a^{t)dW{t)]) (6.79) 

^i{{e\[Hsit),A{t)]dt 

-K{t){{e\ (^A{t)[a^{t)a{t) + a^{t)a{t)]) dt 

+ar(2/t(t)[n(t) + u] + n{t)){{e\A{t)a^{t)a^{t)dt 

+ {{e\ (A{t)[a^{t)dWt + a*{t)dWt]) . (6.80) 

In terms of operators a{t) and a^{t) it becomes 

d{{e\A{t)=z{{e\[Hs{t),A{t)]dt 

-K{t){{e\ {[A{t),a^{t)}a{t) + a^{t)[a{t),A{t)})dt 

-a{2K{t)n{t) + n{t)){{e\[[A{t), a^{t)}, a{t)}dt 

+ m {[A{t), a^{t)}dFt - a[A{t), a{t)}dF^) , (6.81) 

where we used 

{\dWt^{\dFt, {\dWt^r*{\dFl (6.82) 

The stochastic equation of motion of the Stratonovich type for the bra- vector 
state {{9\A{t) is derived similarly in the form 

dmA{t)^i{{9\[Hs{t),A{t)]dt 

-Kit){{9\ (att(t) [a{t) ,A{t)}- aa{t) [a^ (t) , A{t) }) dt 

+ {{9\ {[A{t), att(t)} o dFt - a[A{t), a{t)} o dF}) . (6.83) 

Applying to (6.81) the random force ket-vacuum |) and the ket-vacuum \ 0) of 
the relevant system, one obtains the equation of motion for the expectation 
value of an arbitrary operator A{t) of the relevant system as 



33 



-Km{[m,aHt)}a{t)+aHtMt),A{t)}))) 
-ai2K{t)n{t) + nmHim, «^ W), (6-84) 

Here, ((■ ■ ■)) = (|(^| • • • \ 0)\) means to take both random average and vacuum 
expectation. This is the exact equation of motion for systems with hnear- 
dissipativc coupling to reservoir, which can be also derived by means of the 
Fokker-Planck equation (6.18). Here, we used the property 



{[A(t), aHt)}dFt) = {[A(t), a(t)}dF,^) = 0, (6.85) 

which is the characteristics of the Ito multiplication. Note that equation of 
motion for expectation value of an arbitrary operator A{t) does not depend 
on the parameter A. 



7 Semi-free system with a stationary process 

One possible way to specify a model is to give the Boltzmann equation (3.23). 
For the cases of a semi-free system corresponding to the stationary quantum 
stochastic processes, one needs to make substitutions 

iE<(i) = 2Kn, u;{t) = lo, K{t) = (7.1) 
where n is an average quantum number in equilibrium given by 



n- 



(e-/^-^)"\ (7.2) 



and T is the temperature of environment (here we use the system with the 
Boltzmann constant k-Q — i)- Then, the Boltzmann equation (3.23) becomes 



h{t)^-2K{n{t) -n) . (7.3) 

It describes the system of a damped harmonic oscillator. 

Substituting the Boltzmann equation (7.3) into the semi-free hat-Hamiltonian 
(3.4) with (3.5) and (3.28) or with (3.50) and (3.51), one obtains 



H — uj {a^a — a^a^ -\- 2aTm an) ad + 2aTmna)a^ 
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— ire (1 + 2cm) (a^a + a^a^ — 2iK,n (7-4) 

= ua'^a'' - iKa^A^'^a" + cr(u; + m) (7.5) 

^uj (y^^t - 7^7t) - (7^7* + 7^7*) - 2(7Ti«; (n(t) - n) 7^7*, (7.6) 

where 



^ 1 + 2an -2an 



\2(1 + an) -(l + 2an) 

The Fokker-Planck equation of the model is given by 
d 



(7.7) 



dt 



\0{t)) = -tH\0{t)), (7.8) 



with (7.6). It is solved as (3.53) with the order parameter 

(^|7*7*l 0) = ar (n(0) - n) (l - e-^"*) , (7.9) 

where 7t, 7"^ and their tilde conjugates are defined by (3.47) and (3.40) with 
n{t) being replaced by the solution of (7.3). The expression (3.53) with the 
order parameter (7.9) led us to the notion of a mechanism named the sponta- 
neous creation of dissipation [79,80,81,82,83,84]. 

Introducing a set of new operators 

= (rft^ ^ ^ coUon (d, rd^) , (7.10) 
defined by 

d" = [5^"' , d^" = B'^'a", (7.11) 

with 



B 



^ 1 + an —an 



\ 



-1 1 

the hat-Hamiltonian H can be also written in the form 



(7.12) 



H^u (dU - d^'d) - m {d)d + J^cf) . (7.13) 
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We see that the new operators satisfy the canonical (anti-)commutation rela- 
tion 



= (7.14) 
and that TSC (3.24) for the thermal ket-vacuum \ 0) can be expressed as 

d\0)^T{n{Q) -n)d^O). (7.15) 
It is easy to see from the diagonalized form (7.13) of H that 

d{t) = V-\t) d V{t) = d e-(^'^+'^)*, (7.16) 
d'^(t)^V-\t)d^V(t) = J^e-(")* . (7.17) 

On the other hand, it is easy to see from the normal ordered form (7.6) that H 
satisfies {9\H — 0, since the annihilation and creation operators satisfy (3.48) 
and (3.49). The difference between the operators which diagonahze H and 
the ones which make H in the form of normal product is one of the features 
of NETFD, and shows the point that the formalism is different from usual 
quantum mechanics and quantum field theory. This is manifestations of the 
fact that the hat-Hamiltonian is a time-evolution generator for irreversible 
processes. 

The second law of thermodynamics tells us that for a closed system the entropy 
increment dS of the relevant system should be given by [85] 

dS = dSi + dSe, (7.18) 
dSi>Q, (7.19) 

where dSi is the change of intrinsic entropy of the system, and dSg the change 
due to the heat flow dQ into the system from the thermal reservoir with 
temperature T: 



dS, = ^. (7.20) 



We can check this for the present model [1]. The entropy of the relevant system 
is given by [86] 

S{t) = - {n{t) Inn(t) - a [1 + an{t)] In [1 + an{t)]} , (7.21) 
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whereas the heat change of the system can be identified with 



dQ{t)^ujdn{t) (7.22) 
leading to 



dS,^'^dn{t). (7.23) 

Putting (7.21) and (7.23) into (7.18) for dS and dS^, respectively, we have a 
relation for the entropy production rate [1] 



dSi dS dSe ^ / / s s, n(t)(l + crn) „ 

— ^ = = 2K(n(t) - n) In --^ — ( > 0, 7.24 

dt dt dt ^ ^ ' ' n(l + on{t)) " ' ^ ' 

giving the inequality (7.19). It is easy to check that inequality holds for both 
cases, i.e. n(t) > n and n{f) < n, and that the equality realizes for the thermal 
equilibrium state, n{t) = n, or for the quasi-static process with « — > 0. 



8 Relation to the Monte Carlo wave-function method 

In this section, we will investigate the Fokker-Planck equation (7.8) in order 
to reveal the relation of NETFD to the Monte Carlo wave-function method, 
i.e. the quantum jump simulation [87,88,89,90,91,92,93,94,9] in which evolu- 
tion with a non-hermitian hat-Hamiltonian is described in terms of randomly 
decided quantum jumps followed by the wave-function normalization. 

Let us decompose the hat-Hamiltonian (7.4) as 

^ = + (8.1) 

with 

(a^a - d^d) - m(l + 2an) [a^a + d^d) , (8.2) 

hW = 2iaTK ((1 + an)ad + na^a^) - 2iKn, (8.3) 

and consider an equation: 



^^mt)y^-iH('^\Oo{t)y. (8.4) 
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Note that H^^^ contains cross terms among tilde and non-tilde operators. We 
see that H^^'> and H'^^^ have the properties 

{e\H^^^^-2iK{l + 2an){e\a}a, (8.5) 
{e\H^^^^2iK{l + 2(7n){e\a)a. (8.6) 

Introducing the wave-functions |^(t)) and |^(t)) through the relation 

m))'=m))m)). m 

we have from (8.4) the Schrodinger equations of the form 



and its tilde conjugate, where 



= uja^a -m{l + 2an)a^a. (8.9) 

This procedure is possible because H^^^ does not contain cross terms among 
tilde and non-tilde operators. The Monte Carlo simulations for quantum sys- 
tems are performed for the Schrodinger equation (8.8) [87,88,89,90]. 

The time evolution generated by the hat-Hamiltonian H^^^ does not preserve 
the normalization of the ket-vacuum, i.e. the normalized ket- vacuum \0{t)) 
evolves for the time increment dt as 



{9\Oo{t + dt))' ^ {9\{1 - iH^'^^dt)\0{t)) ^l-dp{t), (8.10) 

with 

dp{t) = i{9\H^^^\0{t))dt = 2n{l + 2an)n{t)dt. (8.11) 

The recipe of the quantum jump simulation is that, for a time increment dt, 

1) when dp{t) < e with a given positive constant e, the normalized ket- 
vacuum evolves as 



\0{t)) ^ \Oo{t + dt))^ 



\Oo{t + dt)y 

1 - dp{t) 
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(8.12) 



2) in the case dp{t) > £, a quantum jump comes in 



Oi{t + dt)) 



iH^^^dt\0{t)) 
dp{t) 



(8.13) 



The time increment dt should be chosen as the condition dp{t) <^ 1 being 
satisfied. 

Averaging the processes \Oo{t)) and \Oi{t)) with respective probabihties 1 — 
dp{t) and dp{t): 



we can obtain the Fokker-Planck equation (7.8). Note that the ket- vacuums 
\Oo{t)) and \Oi{t)) look like satisfying a certain kind of stochastic Liouville 
equation. 



9 Summciry 

The aim of this paper has been to study the system of QSDEs from a physical 
basis. We have formulated everything from the starting point using the method 
of NETFD. In the presented approach, boson and fermion systems are consid- 
ered simultaneously. The obtained results have two fixed parameters: the real 
parameter a specifying different commutation rules for boson and fermion op- 
erators, and the complex parameter r, (2.8), specifying different thermal state 
conditions for boson and fermion systems. Such a combined consideration was 
made possible due to the unification of fermion and boson stochastic calculus 
(Appendix B) , where fermion annihilation and creation processes arc realized 
in a Boson Fock space by means of a simple stochastic integral prescription 
leading to similar multiplication rules for stochastic differentials. 

The dissipation mechanism is considered through the concept of a quantum 
noise, i.e. as a quantum field interacting with the relevant system. In our paper 
we considered two types of interaction with external fields: hermitian (A = 1) 
and non- hermitian (A = 0). With the latter, conservation of the probability 
is satisfied within the relevant system. With the former, information about 
only relevant system is not enough and instead of that we can speak about 
conservation of the probability within the total system: relevant system plus 
environment system. 



0{t + dt)) ^[l-dp{t)]\Oo{t + dt)) + dp{t)\Oi{t + dt)), 



(8.14) 
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As we are concentrated on the stochastic equations, there are two types of 
stochastic calcuhis: Ito and Stratonovich. Correspondingly, equations used one 
or another type of stochastic calculus are classified as QSDE of the Ito or 
Stratonovich types. The Langevin equation of the Stratonovich type (6.65) 
has structure similar to one of the Heisenberg equation of motion for a dy- 
namical quantity in quantum mechanics and quantum field theory. As a re- 
sult of different stochastic multiplication rule, the Langevin equation of the 
Ito type (6.61) contains an extra term proportional to a product of random 
forces dWtdWt- The corresponding Fokker-Planck equation is then obtained 
most easily from the quantum stochastic Liouville equation of the Ito type by 
taking the random average. Though in fermion case the connection with the 
classical Brownian motion is only formal, the Ito/Stratonovich product for- 
mula is the same as in boson case (relations (A. 15) and (A. 16)). The averaged 
equation of motion for a dynamical quantity can be obtained in two ways. 
Prom the Langevin equation by taking both random average and the relevant 
vacuum expectation, or from the Fokker-Planck equation by taking the vac- 
uum expectation of operators corresponding to the dynamical quantity. In our 
study we showed that QSDEs constructed upon hermitian and non-hermitian 
interaction hat-Hamiltonians lead to the same averaged equation of motion 
(6.84) for an arbitrary operator of the relevant system. In the case of station- 
ary semi-free quantum stochastic process, its irreversibility is checked in terms 
of the Boltzmann entropy. We also demonstrated the relationship between the 
presented formulation and the method of quantum jump simulations. 

The approach we followed in this article is rather formal and we are looking 

now for some demonstrative examples of its application for particular prob- 
lems. An interesting result is obtained, for instance, for the model of a con- 
tinuous quantum non-demolition measurement — continuous observation of a 
particle track in the cloud chamber [95,96], and for the system corresponding 
to the quantum Kramers equation [97,98]. More detailed report about them 
will be presented elsewhere. 
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A Ito and Stratonovich calculus 

Definitions of the Ito [24] and Stratonovich [25] multiplications for arbitrary 
stochastic operators Xf and Yf in the Schrodinger representation are given. 
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respectively, by 



XfdYt^Xt{Yt+at-Yt), (A.l) 
dXfYt^{Xt+dt-Xt)Yt, (A.2) 



and 



Xt o dYt = ^{Xt+at + Xt)iYt+at - Yt), (A.3) 
dXt o Yt = {Xt+dt - Xt) \ (Ft+dt + Yt) . (A.4) 

Prom these relations we have the connection formulae between the Ito and 
Stratonovich products in the differential form as 



Xt o dYt = Xt ■ dYt + ^Xt ■ dYt, (A.5) 
dXt o Yt = dXt ■ Yt + ^dXt ■ dYt. (A.6) 

Note that random average of the stochastic multiphcation (A.l) or (A.2) of 
the Ito type is equal to zero. 

Definitions of the Ito and Stratonovich multiplications for stochastic operators 
X{t) and Y{t) in the Heisenberg representation are given in the same form by 



X{t) ■ dY{t) = X{t) [Y{t + dt) - Y{t)] , (A.7) 

dX{t) ■ Y{t) = [X{t + dt) - X{t)] Y{t), (A.8) 

and 

X{t) o dY{t) = - [X{t + dt) + X{t)\ \Y{t + dt) - Y{t)\ , (A.9) 

dX{t) o Y{t) = [X{t + dt) - X{t)] i [Y{t + dt) + Y{t)] , (A.IO) 

where operators X{t) and dX{t) are introduced, respectively, through relations 

X{t) = Vp\t)XtVF{i), (A.ll) 

dX{t) = d{yp^{t)XtVF{t)), (A.12) 

with Vpii) being a stochastic time evolution operator. 
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From (A. 7) to (A. 10), wc have the connection formulae between the Ito and 
Stratonovich products in the differential form as 



X{t) o dY{t) = X{t) ■ dY{t) + lrdX{t) ■ dY{t), (A.13) 
dX{t) o Y{t) = dX{t) ■ Y{t) + \dX{t) ■ dY{t). (A.14) 



Stochastic multiplications (A. 7) to (A. 10) are consistent with corresponding 
types of differential calculus for products of stochastic operators, which for 
the case of the Ito type calculus and the Stratonovich type calculus read, 
respectively, as 

d{X{t)Y{t)\ = dX{t) ■ Y{t) + X{t) ■ dY{t) + dX{t) ■ dY{t), (A.15) 

and 

d[X{t)Y{t)] = dX{t) o Y{t) + X{t) o dY{t). (A.16) 
B Boson and fermion Brownian motion 

Let Fg denotes the boson Fock space (the symmetric Fock space) over the 
Hilbert space IK = L^(R_,_) of square integrable functions, and ht and h\ denote, 
respectively, boson annihilation and creation operators at time t e [0, oo) 
satisfying the canonical commutation relations 

[ht,h\]^5{t-s), [6„6,] = [6t,6t]^o. (B.l) 
The bra- and ket- vacuums (| and |) are defined, respectively, by 

{\h\ = 0, 6,1) = 0. (B.2) 

Note that (| = I)''" since here we are considering the unitary representation of 
bt and 6j. The space F^ is equipped with a total family of exponential vectors 

(e(/)| = (|exp|Jrftr(t)6,|, (B.3) 
\eig)) = expljdtg{t)bl\\), (B.4) 
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whose overlapping is 

(e(/)|e(g)) = exp I Jdt r(t)g(t)| . (B.5) 

Here, f,g G IK. The dense span of exponential vectors is denoted by £. 
Operators bt, bl and exponential vectors are characterized by the relations 

ieiml = (e(/)|r(t), 6,|e(g)) = g(i)|e(g)). (B.6) 
Let us introduce operator C/j defined as 

Ut = 0-<P[o,t] + (7>P(t,oo), (B.7) 

where (j< and (T> are two independent parameters taking values ±1, and P[a,b] 
{a<b) is an operator on !K of multiplication by the indicator function whose 
action reads 

oo b CO 

PlaM J dt g{t) = jdt g{t) = jdt e{t - a)e{b - t)g{t). (B.8) 

a 

Here, 9{t) is the step function specified by 



9{t) = 



1 for t > 0, 
for t < 0. 



(B.9) 



The operator P[a,b] has the following properties: 

Pla,b] = P[a,b]^ P[a,b] = ^[a,6]^[c,d] = ^[c,d]^[a,b] , (B.IO) 

which are easily verified using the definition (B.8). Then, we see that operator 
Uf is unitary, and satisfies 

U^ = I, U} = Ut, UtU, = UMt, (B.ll) 
where / is the identity operator. 

The so-called reflection process Jt = Jt{Ut), t G R+, whose action on £ is 
given by 
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Ji|e(g)) = \e{Utg)) = exp I^Ut J dt' g{t%j |), (B.12) 
inherits properties of the operator Ut, (B.ll), i.e. 

J^^l, jJ^Jt, JtJs^JsJt, (B.13) 
and does not change the vacuum: 

(l^*=(l, ^*l) = l)- (B.14) 
Here, 1 is the unit operator defined in F^. 

Let us now consider new operators 

ht = JA, hi = bUt. (B.15) 
Apparently, they annihilate vacuums 

(|b| = 0, b,|)-0. (B.16) 
The following matrix elements 

(e(/)|[J„6,]_.|e(g)) = (e(/)|{l 

-a[a> + (a< - a>)^(i - s)]}g{s)Jt\e{g)), (B.17) 
(e(/)|[b„bt]_,|e(g)) = (e(/)|5(t-.)|e(g)) 

+ (e(/)| J,J,(a>a< - a)/* (s)g(t)|e(g)), (B.18) 

are valid for f,gEJi. Then the requirement of equal-time (anti-)commutati- 
vity between Jf and bt 

[Jt,bt]-a^O (B.19) 

gives 

l-aa< = 0, (B.20) 
while the requirement of canonical (anti-)commutation relation 

[ht,hlU = S{t-s) (B.21) 
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leads to 



(T>(7< — a — 0. 



(B.22) 



All those conditions are satisfied when — a and (7> = +1. Then the 
operator Ut turns out to be 



Ut = crP[o,t] + P{t,oo)- 



(B.23) 



Note that for a boson system, i.e. a — 1, Ut ~ I and the operators and bj 
reduce, respectively, to bt and 6j. 

We see that the generalized quantum Brownian motion, defined by 



Qt = J dt' b*/, B| = j dt' bj,. 



(B.24) 



with Bo = 0, Bq = 0, satisfies 



[Bi,Bt]_, = min(i,s). 



(B.25) 



The case a = 1 represents the boson Brownian motion [49,55], whereas the 
case cr = — 1 the fermion Brownian motion [62]. Their increments 



dBt = Bt- 



-dt 



Bt = btdt, 



dBl = Bl,, - Bl = hldt, 



annihilate the vacuum, i.e. 



[B.26) 
(B.27) 



{\dBl = 0, dBt\)=0, 
and their matrix elements read 



(B.28) 



{e{f)\dBt\e{g) 
{e{f)\dBl\eig) 
{e{f)\dBtdBt\e{g) 
{e{f)\dBldBt\e{g) 
{e{f)\dBtdBl\e{g) 



= {e{f)\Jtg{t)dt\e{g)), 
= {e{f)\r{t)dtJt\e{g)), 
= 0, 
= 0, 

= dt{e{fMg)). 



(B.29) 
(B.30) 
(B.31) 
(B.32) 
(B.33) 
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Here we neglected terms of the higher order than dt. The latter equations are 
summarized in the following table of multiplication rules for increments dBt 
and dBl: 



dBt 

dBl 
dt 



dBt dBl dt 



dt 





(B.34) 



Now we consider a tensor product space f = r° f°. Its vacuum states |)) 
and exponential vectors |e(/, g))) are defined through the "principle of corre- 
spondence" [3] 



I))— (B.35) 
|e(/,g)))"|e(/))(e(g)|. (B.36) 

Annihilation and creation operators acting on F are defined through 



bt\e{f,g)))^bt\e{me{g)\, (B.37) 

6l|e(/,g)))^6l|e(/))(e(g)|, (B.38) 

bt\e{f,g)))^\e{me{g)\bl (B.39) 

bl\e{f,g)))^\e{me{g)\bt, (B.40) 

and similarly for Jt and Jt, i.e. 

J,|e(/,g)))«J,|e(/))(e(g)|, (B.41) 

Jt\e{f,g)))^\e{f)){e{g)\Jt. (B.42) 

Algebra of commutation relations between these operators reads 

[bt,bl]-[bt,bl]^S{t-s), (B.43) 

[btM-[Kbl]-0, (B.44) 

[Jt,6,] = [Jt,6,]=0, (B.45) 

[Jt,btU^[Jt,btU^Q- (B.46) 



Let us now consider new operators defined by 
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ht^Jth, bl^bUt, (B.47) 
bt^rJA, bl^rbU,, (B.48) 

where f is an operator satisfying the following (anti-) commutation relations 



[f, Jt] = [r, .It] = 0, (B.49) 

[f,btU = [f,blU = 0, (B.50) 

[f,btU = [f,blU^O, (B.51) 

and the condition 

= a. (B.52) 
Operators bt, bj and their tilde conjugates annihilate vacuums 

((|bl = ((|bl = 0, b,|)) = b,|)) = 0, (B.53) 
and satisfy canonical (anti-)commutation relations 

[bt,blU = KblU = S{t-s), (B.54) 

[bt,b,]_, = [bt,bl]_, = 0. (B.55) 

Since (bt)t and (bl)~ are calculated as 

{ht)''^blJtr^ ^ar^bUt, (B.56) 

{blr^fibUtr^fblJt, (B.57) 

the commutativity of tilde conjugation and hermitian conjugation for opera- 
tors (B.47) and (B.48) imphes 



= af. (B.58) 

In order to fulfill the requirement that double tilde conjugation applied to 
operators b^'s leaves them unchanged one needs to put 



(f)- = ft, (B.59) 



smce 



(bt)~ = f(f)~JA = f(f)~bt. (B.60) 
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Because of 



bt|)) = &t(r|)))=0, (B.61) 

one can conclude that f|)) oc |)). The proportionahty factor is a phase factor 
since the norm of f|)) is unity: 



((|ftf|)) = ((|aff|)) = a' = l. (B.62) 
Hence one can write 

f|)) = e^'^/^|)). (B.63) 
Multiplying both sides by f, one has 

f'\))=e^^\)) = a\% (B.64) 
which gives e^'^ — a, or 

f|)) = v^|)). (B.65) 

Thermal degree of freedom can be introduced by the Bogoliubov transforma- 
tion in f . For this purpose we require that the expectation value of bjb^ should 
be 

(bib,) = n5{t - s) (B.66) 

with n e R+, where (■ ■ ■) = (| • ■ ■ |) indicates the expectation with respect to 
tilde invariant thermal vacuums (| and |). The requirement (B.66) is consistent 
with TSC for states (I and 1) such that 



|bl = r*(|b„ b,|) = ^^bl|). (B.67) 



Let us introduce annihilation and creation operators 

Q = [1 + an]ht - arnhl, (B.68) 

cf^hl-arht, (B.69) 

and their tilde conjugates. Prom the TSC (B.67) one has 



48 



{\cf^{\ct^O, q|) = q|) = 0. (B.70) 
With the thermal doublet notations 

= (bl, -Tht) , br = coUon (bi, rbl) , (B.71) 

and 

cf = (c*, -tq) , cr = coUon (q, rcf) , (B.72) 

(B.68), (B.69) and their tilde conjugates can be written in form of the Bogoli- 
ubov transformation 

c^ = B^''hl -c'^ = hnB-T, (B.73) 

with (7.12). This new operators satisfy the canonical (anti-)commutation re- 
lations 



[ct,ctU = 6{t-s). (B.74) 

In the following, we will use the representation space constructed on vacuums 
(I and I). Note that (| ^ \)\ i.e. it is not a unitary representation. Let 
denotes the Fock space spanned by the basic bra- and ket-vectors introduced 
by a cyclic operations of q, q on the thermal bra- vacuum (|, and of cf , cf on 
the thermal ket- vacuum |). Quantum Brownian motion at finite temperature 
is defined in the Fock space by operators 



J ds bl Bl = Jds bl (B.75) 







with Bq = and Bq = 0, where jj stands for null or dagger. The explicit repre- 
sentation of processes Bf and Bf can be performed in terms of the Bogoliubov 
transformation. The couple B^ and bJ, for example, is calculated as 



t 

Bt 



t 



■ j ds {Cs + arncf) = Q arnCf, (B.76) 



t 

Bl = Jds{[l + an]ct + tQ = [1 + an]Cf + tQ, (B.77) 



where we defined new operators 
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I 



ds c«, C« 



z 

I 



ds cl 



(B.78) 



with Cg = and Cq = 0, and jj standing for null or the Venus-mark. Since 
matrix elements of dCl and dCl in thermal space f read 



{dCt) = {dQ) = {dCf) = {dCf)^0, 
{dCfdCt) = (dCfdCt) = 0, 
{dCtdCf) = {dCtdCf) = dt, 



(B.79) 
(B.80) 
(B.81) 



calculation of moments of quantum Brownian motion in the thermal space 
can be performed, for instance, as 



(dBtdSl) = ((dCt + aTndCf){[l + an]dCf + rdCtj) 

= [1 + an]{dCtdCf) = [1 + an]dt. 



(B.82) 



Repeating this for other pair products of rfBf, rfBf and dt, multiplication rules 
for these increments can be summarized in the following table: 



dBt 

dBl 

dBt 



dBl 
dt 



dBf 



dBl 



dBf 



dB 



t 



dt 



[1 + an]dt Tfidt 
ndt T[l + an\dt 

arndt [1 + (Tn]dt 

aT[l + an]dt ndt 




(B.83) 



C Treatment of fermions in Thermo Field Dynamics 



We arc deciding the double tilde conjugation rule and the thermal state con- 
ditions for fermions [99,100] by considering the system consisting of a vector 
field and Faddeev- Popov ghosts [101]. 

In the case of pure Abelian gauge field within the Feynman gauge, the system 
is specified by the Hamiltonian -f^vf+gh = -^vf + -f^gh defined on the total state 
vector space V = Vvf Vgh- and ifgh are, respectively, Hamiltonians for 
the vector field and ghosts defined on the vector field sector Vyf and the ghost 
sector Vgh given by 
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(C.l) 



— * — * 

with e{k) — \k\ being the energy spectrum and g^'^ = diag(l, —1, —1, —1), and 

by 



i/gh ^-i j d^k e{k) [c\k)c{k) - c\k)c{k) 



(C.2) 



, — * — * 

Here, aJ((A;) and a^(A;) are, respectively, creation and annihilation operators of 

— * 

the gauge field of the mode k satisfying the canonical commutation relations 



K(^), alio)] = -g^J^ik - q), 



(C.3) 



while c'lk) and c{k) [c'lk) and c{k)] are, respectively, creation and annihila- 
tion operators of ghosts [anti-ghosts] satisfying the following canonical anti- 
commutation relations: 



(C.4) 



Other combinations of ghost/anti-ghost operators anti-commute with each 
other. The BRS charge - generator of the BRS transformation [102], and the 
ghost charge [103,104] acting on the total space V are, respectively, given by 



Qb^- J d^k k^ [a^{k)c\k) + al{k)c{k) 



Qc 



c\k)c{k) + c\k)c{k) 



(C.5) 
(C.6) 



which satisfy 



[iQc, Qb] = Qb- 



(C.7) 



Let us introduce a set of new operators {a(a){k^)\(^ = +, —, L, S} through the 
relation 



af,{k) = a(a){k)el^\k), 
where e^^^^k) are polarization vectors defined by 



(Ci 



6W(fc) = (0,e±), 



(C.9) 
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ef{k)^-ik,m''^-i{\k\-k)/2\k\^ 



(C.IO) 
(C.ll) 



with e± satisfying etj. ■ = 0, ■ Cip = and ■ e± = 1. The polarization 
vectors e^^\k) correspond, respectively, to the transverse modes with helicity 
±1, while e[i^''(A;) and £^^\k) indicate, respectively, the longitudinal mode and 
the scalar mode. With the definition (C.IO) and (C.ll), we see that e^^>{k) ■ 
e(L),/^(A^) = e(f)*(fc) • e(S)''^(A^) = 0, and )*(fc) • e^^^'^'ik) = 1. Introducing a 
"metric" 



V 



-1 o^ 

0-100 
1 
1 Oy 



(C.12) 



we can define "contravariant" polarization vectors through e^^^ (k) — g(^^^-^e^'^''^{k), 
and see that 



ei:^{k)-e^'^^^''ik) = r, 
e^^>{k) ■e^^^'^'{k)^g^^^\ 



(C.13) 
(C.14) 



The commutation relations (C.3) being rewritten in terms of operators a(^^){k) 
and ci((^)(/c) become 



[a^a){k),al^^{q)] = -g(^r)S{k - q). 



(C.15) 



Also the Hamiltonian for the vector field and generator of the BRS transfor- 
mation read 



^vf = -/ d'ke{k)^''^^al^{k)a^,){k), 
Qb^-i j d^k [a\s-^{k)c{k) - a(s){k)c\k) 



(C.16) 
(C.17) 



In the local covariant operator formalism [103,104] of gauge theories, the space 
V of state vectors has inevitably an indefinite metric as can be seen by (C.15) 
with (C.12). The physical subspace Vphys of V, defined by [104] 



(C.18) 



52 



can be shown to have a positive semi-definite metric [104]. Dividing Vphys by 
its subspace Vq consisting of normless states, we have, as a quotient space, 
the physical Hilbert space i^phys (= Vphys/Vo) with positive definite metric 
in which the probabihstic interpretation of quantum theory works, -ffphys is 
isomorphic to the Hilbert space ^Kphys spanned by the Fock states created by 
the cyclic operation of a|j_^(^) on a certain vacuum. The space spanned by 
the Fock states created by a|s)(^) is classified in Vo, while the space spanned 

by the Fock states created by aj^) (k) is classified in a space complemented to 
Vphys- This reflects the fact that physical modes for photons are two transverse 
modes only. 

As it is sufficient to pay attention to one mode in the following manipulation, 
we will pick up a mode, say k, from each type of particles, and drop the index 
k, for simplicity. Let us span the state vector space V by means of a set of the 
bases {|{n(o-)}) ■ |nc, ng)} whose elements, being the bases of the vector field 
sector Vvf and the ghost sector Vgh, respectively, are defined by 



n ^ I {0})' (c-19) 

o-=±,L,S V^(a)' 



|ne,n,) = (ct)"^(ct)"'|0,0). (C.20) 
They constitute the eigenstates of i^yf, -f^gh and iQc- 

ff.f|{n(.)}) = E.f({^w})|{^w}), (C.21) 
^ghl^^c^^c) =^gh(nc,ng)|nc,ng), (C.22) 
iQc\nc,nc)=Ngh{nc,nc)\nc,nc), (C.23) 

with Eyf{{n(„)}) = £T,a'n{cy), Egt,{nc,nc) = £ (ric + ng) and iVgh(nc,ng) = 
{ric — nc), where n(o-) are non-negative integers, and ric and ng take values 
of or 1. We will denote the basis vectors {|{n(o-)}) ■ |nc,ng)} by {|n)} for 
brevity. Then, the metric tensor of V is 

Vn,m = {n\m). (C.24) 



At the finite temperature, the statistical average of an observable quantity A, 
satisfying 



[Qb,A]=Q, (C.25) 
is given by [105] 



53 



(A) =Tr ApP(°) = Tr Ape''^'^ 



(C.26) 



with P*^°) being a projection operator onto IKphys and p = e~^'^^f+8'i be- 
ing the statistical operator acting on V with the partition function Z — 
g-^i?vf+gh+7rQc. ^]^g trace operation is taken in the space V. Here, for the 
second equahty in (C.26), we used the BRS-invariance of the statistical oper- 
ator 

[Qb,p]^0. (C.27) 

Let us express the statistical average (C.26) as the vacuum expectation in the 
doubled state space (thermal space) V = V (8) V which is introduced as follows. 
If A is an operator on V so that 

A = J2Anm\n){ml (C.28) 

n,m 

the corresponding vector 1^4) in V is obtained as 

\A) = A^^\n,m), (C.29) 

where {\n, fh) = \{n(^a-)}, {^(o-)}) ■ l^o 't-c; t^c, "^c)} is the set of the bases span- 
ning V and defined through the "principle of correspondence" [3] : 

\{n(a)}, {rri{a)}) < — > |{^M})({^(a)}|, (C.30) 
|nc, Uc, rhc, rhc) < — > |nc, ric) (mc, mc| . (C.31) 

The inner product in V is given by 

{A\B) = Tr A'^B. (C.32) 
Annihilation and creation operators acting on V are defined through 



. ^ \ \{'n(a)}){{m(a)}\, (C.33) 



\{nia)},{^{a)}) ^ \{n{a)}){{m^a)}\ 



{CM) 
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for vector field, and through 



\nc,nc, rhc,rhc) 



|nc,nc)(mc,mc|, 



(C.35) 



|nc,nc,mc,mc) 



|nc, nc)(mc,mc|, 



(C.36) 



\nc,nc,mc,mc) 



1)"+ \nc,nc){mc,mc 



(C.37) 



\nc,nc,rhc,rhc) (-l)"|nc, nc)(mc, ma 



(C.38) 



for ghosts, where u = Ngh{nc,nc) — A^gh(mc, mg). Also bases \{'n(a)} , {fh(a)}) 
and \nc, ric, rhc, rhc) are generated from the vacuums |{0}, {0}) and |0, 0, 0, 0), 
respectively, as 



|{0},{0}), (C.39) 



\n,,n„m,,m,) = (-l)-^'(c^)"1ct)"^(ct)™^(gt)'«^|0, 0, 0, 0), (C.40) 
where v — mcrric- For the total vacuum, we will use a collective designation 

|o,6). 

Let us introduce thermal vacuums (^| and |f £ V such that 



{A) = {e\A\om- 

We require them to satisfy 



(C.41) 



(C.42) 



and 
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(^|4 = 0, Q,\0{P))^0, (C.43) 

where Qj^ and Qc are the generator of the BRS transformation and the ghost 
hat-charge, respectively, in V [99] 



Qb = Qb - Qb = -i [a\sf - a(s)C^ + aJs)C - a(s)C^] , (C.44) 
Qc = Qc- Qc = c^c + c^c- c^i-c^c. (C.45) 

To satisfy (C.42), we need a trick. Namely, by rewriting (C.26) as 

(A) = Tr Mpe'^^-rS (C.46) 

we introduce an operator 9 with the basic requirement that its inverse exists. 
Then we settle the correspondence 

(^1^^, \0{(3)) ^ pe'^^^e-K (C.47) 
It gives 

(^IQ^ = (^1 (c^c + c^c - c^i - c^c) 

^ e (c^c + c'^c) - (c^c + c^c) e 

= [9,Qc]. (C.48) 

With the requirement given by the first equality in (C.43), the expression in 
(C.48) is equal to zero and tells us that 9 and Qc commute with each other, 
i.e. [9, Qc] = 0. Then 



Qc\0{/3)) = [ck + ck-c^^-^^c) \ 0{l3)) 



+ c^c)p e^^^ 9-^-p e^^^ ^"^ (c^c + c^c) 



c'c 



- [gc,p]e"«^r\ (C.49) 

and the second equahty in (C.43) is automatically satisfied as far as [Q^ p] — 0- 

Based upon [9,Qc] = and existence of 9^^, let us try the foUowing form of 
9: 

where, 0i, 02 and 03 are real numbers which should be determined. Then, 
taking into account the first correspondence in (C.47), the definition (C.44) 
and (C.50), the calculation of the first equality in (C.42) goes as 
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{9\Q'b^ -i{9\ (a{s)C - a(s)C^ + ajg^c - a(s)c'^) 

^ -i (a|g)C - a(s)C^) + (a(s)C^ + a(s)c) 6'} 
= -ie {a{s)C - a(s)C^ + e-''^i-^'^3(-ic-tc) ^^^^^t 

If we take 0i = 0, 02 = ti" and 03 = 0, we have (6'|(5b = with the choice 

e = exp{m{-ic^c)}. (C.52) 
This structure for 9 allows us to calculate the second equahty in (C.42) as 

Qb\0{I3)) = ~t (a{s)C-a(s)c"^ + a{g)C-a(s)c'f) \0{P)) 

= [gB,p]e"«^r\ (C.53) 

where we also used [e'^'5=, Qb]+ — which is obtained from (C.7). Taking into 
account the BRS-invariance of the statistical operator p (C.27), the expres- 
sion (C.53) is equal to zero, and both requirements (C.42) are fulfilled with 
the choice (C.52). We see that introduction of factors 9 and 9~^ indeed is nec- 
essary to satisfy the BRS-invariance of the thermal vacuums (^| and \0{/3)). 
Expression of the unit operator in V 

l = El^Km(^i|, (C.54) 

n,m 

and correspondences (C.47) with (C.52) enable us to see the structure of ther- 
mal vacuums as 



{0\ = ^{Vnlnn,rh\0 



E iVnl^,m^^/{{n{a)},{^{a)}\ 

"(<t).™(<t) 

E (^rne,ne-),(me,mc-))*K, U,, rhc, mc\9 



X 



nc,nc 



({0}, {0}| cxp {-g(--)5(,)a(,)} (0, 0, 0, 0|[1 + icc\[l + ilc] 
(0, 0| exp jicc -|- ice — g^'^^''a(o.)a(r)} , 



(C.55) 



and 
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|0(/3)) = Z-^exp{-e-^^g(^^)a}^)a[^) -ie-^^ (c'd^ + d^c')} |0,0). (C.56) 
It may be instructive to note that considering 



(C.57) 



with Qb = Qb + Qb, instead of (C.42) with (C.44), leads to the choice 9 = 
exp{iTr{ic^c)}. 

After determination of parameters 0i the thermal state conditions with the 
ghost operators can be derived through the following steps. First, for the bra- 
vacuum we see 







(A 















—c 



(C.58) 



-ct 



therefore 



(^1 



(A 













(C.59) 



Similarly, taking into account structures of the statistical operator p, the ghost 
charge Qc, (C.6), and 9, (C.52), for the ket- vacuum we have 



/c^ 



/ e-/5^ct 



-e-^^gt 



\0{P)) 



\0{f3)) 



p e^^<= 9-^ 



■p e 



TtQc Z]— 1 



- e-^^ c 



which gives 

/ 
V 



0{(3))=e 



-f3e 



( ct 



-gt 



(C.60) 



(C.61) 



58 



The double tilde conjugation rule must be defined so that it leaves thermal 
vacuums unchanged. To this end, we put 




and determine parameters ^ and ^ so that (^^|'" 
Taking the tilde conjugation of (C.55) we have 



(C.62) 
(C.63) 

(^1 and|0(/3))- = |0(/3)). 



(6i|~ = (0, 0|~ exp{-i^cc - i^cc - g('"^'a(^)a(^)} 



= (0, 0| exp{i^cc + i^cc - S-''^^ai^a)a{T)}-, 

where we assumed (0, 0|~ = (0, 0|. The requirement {Q\ 
^ = 1, which leads to 



(C.64) 

= (^1 gives ^ = 1 and 













V) 



(C.65) 



As a consequence, we obtain the tilde-invariance of the thermal ket-vacuum 
too: 



(z-^exp{-e-^^g(-)aj,)a|,) + ctct)}|0, 0))^ 

= Z-^ exp{- e-^^ g("")5[,)a;,) + i ^'^'{^c^ + 0) 
-W)). (C.66) 



A similar line of reasoning can be used to derive the tilde conjugation rule and 
the thermal state conditions for a system consisting of physical fermions. Let 
us consider the system specified by the Hamiltonian 



d^keik) a^{k)a{k) 



(C.67) 



— * J. ~* 

with a{k) and a'{k) being, respectively, fermion annihilation and creation 
operators satisfying the canonical anti-commutation relation 



[a{k),a\q)]+ = 6{k - q). 



(C.68) 
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As in previous consideration, in the following manipulation wc will pay atten- 

— * — * 

tion to one mode, say k, and drop the index k for simplicity. In that case the 
bases of the state vector space Vf will be denoted as |0) and |1) defined by 
a|0) = and |1) =at|0). 

At the finite temperature, the statistical average of an observable quantity A 
is given by 



(A) = Tr Api 
with pf being the statistical operator of the system 



(C.69) 



p, = Zi^e-^"^ = [|0)(0| +e-^^ |1)(1| 
where Zf is the partition function. 



(C.70) 



Within the doubled state space Vf = Vf Vf, the statistical average (C.69) 
can be expressed in terms of the vacuum expectation with respect to the ther- 
mal bra- and ket-vacuums for physical fermions. The bases of Vf are defined 
through the principle of correspondence: 



n, m) 



|n)(m|. 



(C.71) 



where n and m take values of or 1. Annihilation and creation operators 
acting on Vf are defined through 







(a\ 


(A 








n, fh) < — > 


[a)) 





\n){m\, 

a\n, fh)i — ^(-l)"-™+V)(m|a^ 
a^|n, fh) < — (-l)"~'"|n)(m|a. 

The bases \n, rh) are generated from the vacuum |0, 0): 

\n,m) = (a^)"(a^)'"|0,0). 



(C.72) 

(C.73) 
(C.74) 



(C.75) 



This time we don't have peculiar symmetries for thermal vacuums like in the 
case of gauge theories with the BRS symmetry. However, in the derivation of 
thermal vacuums, let us use a trick similar to the one in (C.46) 



(^) =Tr e^-^^^^Ap e-^"^"^" 



{e\A\om, 



(C.76) 
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where is a real number which should be decided. We settle the correspon- 
dence for (^1 and \ 0{P)) as 



(^1 



< — > e 



\0{P))^p,e 



(C.77) 



They are normalized, {9\0{(3)) = 1, and generated from |0,0) as 



1 + aa 
l + e-'-^e-^^a^a^ 



|0,0). 



(C.78) 
(C.79) 



A requirement of the tilde invariance for the thermal bra- vacuum 



(^1 



(C.80) 



determines the tilde conjugation rule for physical fermion operators up to the 
phase factor 



(a)~ = a, (a)~ = -e^^-^a. 



(C.81) 



We have seen in (C.65) that the ghost operators, which are fermion operators, 
were unchanged under the double tilde conjugation. Let us adopt the same 
rule for physical fermion operators, i.e. put = 7r/2 to obtain 

(a)~ = a, (ay = a. (C.82) 
With this choice of 0, the thermal vacuums for physical fermions read 



(^| = (0,0| [1 + iaa], 
\0{(3))^Zf^ [l-ie-^'a^a^] |0,0) 

and satisfy the following thermal state conditions 



(C.83) 
(C.84) 



|at = -i{e\a, a\0{(3)) = ie-'^'a^O{p)). 



(C.85) 



D Correlation of random force operators 



The random force operators are of the Wiener process whose first and second 
moments are given by real c-numbers: 



61 



{dFt) = {dFl) = 0, (D.l) 

{dFtdFt)^{dF}dF})^Q, (D.2) 

{dFtdF}) — a real c-number, (D.3) 

(dF^dFt) = a real c-number, (D.4) 

where (■ ■ ■) = (| ■ ■ ■ |) represents the random average referring to the random 
force operators dF^. From (D.l), (D.2) and TSC (6.4) we have for operators 
(6.14) and (6.15) 

{dWt) = {dWt) = {dWD = {dWt^) = 0, (D.5) 
{dWtdWs) = {dWtdWs) = 0, (D.6) 

while from (6.16) it follows 

{dWldWs) = {dWldWs) = 0, (D.7) 
{dW^dW^} = {dW^dW^} = 0, (D.8) 



and their tilde conjugates. Using (D.6) to (D.8) the exphcit structure of 
dMfdMt in (6.23) is written as 

dMtdMt = -2adWtdWta^a* + A (dWtdW^^a^a + dWtdWt^a^a) (D.9) 

in a "weak sense" ^ . We demand that the Stratonovich type time evolution 
generator should not contain a diffusion term, i.e. the term proportional to 
a^a^. Then the correlation {dWtdWt) is determined to be 

dWtdWt = T {2K(t)[n(t) +r]]+ h(t)} dt = (dWtdWt) (D.IO) 

so that 77d in (6.23) is cancelled by the first term in the r.h.s. of (D.9). Here, 
the first equality in (D.IO) should be understood in a weak sense as well. 
Expression (D.IO) is compatible with the assumption that the process is white. 
Let us put the subscript F to S^(t) in the Boltzmann equation in order to 
remember that it is due to the interaction with the random force dF^. 

h{t)^-2K{t)n{t)+iJ:^{t). (D.ll) 
Making use of two previous equations, we have 

^ In the case of classical systems it corresponds to the stochastic convergence. 
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i^f{t)dt = 2K{t)n{t)dt + h{t)dt 

= -2K{t)r]dt + (TT{dWtdWt) 

= -2K{t)r]dt + {dF^dFt) + i/ [{dFtdF^t) - (7{dF}dFt)\ , (D.12) 
where (6.14) has been used, and has been erased with the help of (5.32). 

We can assume that the quantity r) may depend on u, i.e. 7] ~ ^(i^), and that 
the physical quantities K,{t), S^(t), (dFfdFt), and {dFtdP}) may not depend 
on u. Then, differentiating equation (D.12) with respect to i/, one has 

= -2K{t)^dt + {dFtdFl) - a{dF^dFt). (D.13) 
This leads to 

drj 
dv 

which is solved as 



m (D.14) 



r]^k{t)u + l{t), (D.15) 

where k{t) and l{t) are real numbers independent of z/. With this solution one 
has 



{dFtdP}) - o{dF}dFt) = 2K{t)k{t)dt, (D.16) 

and 

iJ:f{t)dt^-2K{t)l(t)dt + (dF^dFt), (D.17) 
which leads to 

{dF}dFt) = {2K{t)[l{t) + n{t)] + h{t)} dt, (D.18) 
where we have used (D.ll). The substitution of (D.18) into (D.16) gives us 

{dFtdp}) = {2K{t)[k{t) + al{t) + an{t)] + ah{t)} dt. (D.19) 



For the case of stationary quantum stochastic process, the Boltzmann equation 
(D.ll) reduces to 
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n{t) = -2K[n{t) - n], (D.20) 

where n is the average quantum number in equihbrium. Therefore, (D.18) and 
(D.19) reduce, respectively, to 



{dF}dFt) = 2K[n + l{t)]dt, (D.21) 
{dFtdF}) = 2K[k{t) + al{t) + an]dt. (D.22) 

Since in the white noise assumption the Boltzmann equation (D.20) is com- 
patible with the stationary process specified by [41] 

{dF}dFt) = 2Kndt, (D.23) 
{dFtdF}) = 2k[1 + an]dt, (D.24) 

one concludes now that 

l{t) = 0, k{t) = 1, (D.25) 
which leads to 

7] = u, ^ = fi. (D.26) 

Note that the result (D.24) can be obtained using the Bogoliubov transforma- 
tion as it is described in Appendix B. 

Substituting (D.25) into (D.18) and (D.19), one obtains 

{dF}dFt) = [2K{t)n{t) + h{t)] dt, (D.27) 
(dFtdFt^) = {2K{t)[l + an{t)] + ah{t)} dt, (D.28) 

which leads to 

{dWtdWt^) = {dFt, dF}) - a{dFt\ dFt) = 2K{t)dt. (D.29) 

Assembling (D.9), (D.IO), (D.26) and (D.29) one obtains expression (6.25). 
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